Charge and momentum transfer in supercooled melts: 
Why should their relaxation times differ? 
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The steady state values of the viscosity and the intrinsic ionic-conductivity of quenched melts are 
computed, in terms of independently measurable quantities. The frequency dependence of the ac di- 
electric response is estimated. The discrepancy between the corresponding characteristic relaxation 
times is only apparent; it does not imply distinct mechanisms, but stems from the intrinsic barrier 
distribution for a-relaxation in supercooled fluids and glasses. This type of intrinsic "decoupling" 
is argued not to exceed four orders in magnitude, for known glassformers. We explain the origin 
of the discrepancy between the stretching exponent j3, as extracted from e(w) and the dielectric 
modulus data. The actual width of the barrier distribution always grows with lowering the temper- 
ature. The contrary is an artifact of the large contribution of the dc-conductivity component to the 
modulus data. The methodology allows one to single out other contributions to the conductivity, 
as in "superionic" liquids or when charge carriers are delocalized, implying that in those systems, 
charge transfer does not require structural reconfiguration. 



I. INTRODUCTION 

Molecular motions in deeply supercooled melts and 
glasses are cooperative so that transporting a single 
molecule requires concurrent rearrangement of up to sev- 
eral hundreds of surrounding molecules. Such high de- 
gree of cooperativity results in high barriers even for the 
smallest scale molecular translations. These high barri- 
ers underly the slow, activated dynamics in deeply super- 
cooled melts and the emergence of a mechanically stable 
aperiodic lattice, if a melt is quenched sufficiently rapidly. 
The Random First Order Transition (RFOT) method- 
ology, developed by Wolynes and coworkers, provides a 
constructive microscopic picture of the structural rear- 
rangements in supercooled melts and quenched glasses. 
The RFOT has quantitatively explained or predicted the 
signature phenomena accompanying the glass transition, 
including the connection between the thermodynamic 
and kinetic anomalies [3, 0, H , the length scale of the 
cooperative rearrangements [2|, deviations from Stokes- 
Einstein hydrodynamics aging [f|, the low tempera- 
ture anomalies |,0,|, and more. (See 0| for a recent 
review.) 

Perhaps the most dramatic experimental signature of 
the glass transition is the rapid super- Arrhenius growth 
of the relaxation times with lowering the temperature, 
from about a picosecond, near the melting point T m , to 
as long as hours, at the glass transition temperature T g . 
These relaxation times are deduced via several distinct 
experimental methodologies and all display an extraordi- 
narily broad dynamical range. Nevertheless, making de- 
tailed comparisons between those distinct methodologies 
has required additional phenomenological assumptions. 
Mysteriously, these comparisons show a significant de- 
gree of mismatch, sometimes by several orders of magni- 
tude. For example, the phenomenological "conductivity 
relaxation time" t ct , is consistently shorter than the 
mechanical relaxation time r s , especially at lower temper- 
atures. The apparent time scale separation varies wildly 



from system to system: For instance in molten nitrates, it 
is about four orders of magnitude at T g [TT[ | , while in sil- 
ver containing superionic melts, the r s /r CT ratio becomes 
as large as 10 11 i.e. almost as much as the whole 
dynamical range accessible to the melt! This disparity 
suggested that the mechanical relaxation and the electri- 
cal conductivity in these systems were in fact due to dis- 
tinct mechanisms: At higher temperatures, the time scale 
separation is small so that the two processes strongly af- 
fect each other, or "mix" , while at lower temperatures, 
the processes become increasingly "decoupled" [l3j]. At 
such low temperatures, the mechanical relaxation occurs 
via the aforementioned, activated concerted events, also 
called the primary, or a-relaxation. Other processes that 
seem to decouple from the mechanical relaxation include 
nuclear spin relaxation, rotational diffusion, and the dif- 
fusion of small probes. (For reviews, see [HI, Q EH ) 

Here we focus on two specific transport phenomena: 
low-frequency momentum transfer, i.e. the viscous re- 
sponse, and the ionic conduction in supercooled melts. 
Notwithstanding the complications needed to analyze the 
electrical modulus data [jj, [13, GJ, [H, [H HH E|| , the 
mismatch between the typical relaxation times, corre- 
sponding to the two types of transport, is clearly present. 
Furthermore, in the case of superionic compounds, one 
may show (see below) that conduction occurs without 
distorting the liquid's structure beyond the typical vibra- 
tional displacements. This is much less obvious for com- 
pounds where the ionic motions are "decoupled" from 
the bulk structural relaxation by four orders of magni- 
tude or less, the latter dynamic range comparable to the 
breadth of the a-peaks in dielectric dispersion in insu- 
lating melts near T g (see e.g. [Hj]). Accounting for the 
distribution width is essential here because the viscosity 
and conductivity are distinct, in fact exactly reciprocal 
types of response: In momentum transfer, the velocity 
gradient is the source, and the passed-on rate-of-force is 
the response. In charge transfer, the force on the ion is 
the source, while the arising velocity field is the response. 
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Consistent with this general notion, in computing the 
viscosity, we will average the relaxation time, with re- 
spect to local inhomogeneities, while the intrinsic ionic 
conductivity will be determined by the average rate of a- 
relaxation. Because of the mentioned, extremely broad 
distribution of structural relaxation times r, the quantity 
(t) may reach several orders of magnitude, and so 

an apparent decoupling is indeed expected; no additional 
microscopic mechanisms need to be invoked. 

The microscopic calculation and comparison, of the 
viscosity and the ionic conductivity, are thus the main 
focus of this article. The two quantities are computed, 
in terms of the barrier distribution and other measurable 
material properties, in Sections[TT]and|TTT]respectively. To 
perform comparisons with experiment and assess the up- 
per limit on the "inherent decoupling" between the two 
phenomena, we will discuss the barrier distribution in 
some detail, in Section [TV] We will find that indeed, the 
degree of decoupling should increase with the width of the 
barrier distribution, and hence at lower temperatures, as 
demonstrated by the RFOT methodology [|[ . We will as- 
sess the deviation between a-relaxation times, as deduced 
from viscosity, ionic conductance, and the maximum in 
e"(uj). Further, we will exemplify potential ambiguities in 
using the dielectric modulus formalism in estimating the 
relaxation time distribution. The latter techique has sug- 
gested that for some substances, the distribution width in 
fact decreases with lowering the temperature, in conflict 
with the present results and the correlation between the 
stretching exponent [3 and temperature, predicted earlier 
by the RFOT theory || . We will see that the conflict is 
artificial and results from the large contribution of the ac- 
component to the modulus data, consistent with earlier, 
phenomenological arguments [lT], [HI . 



II. VISCOSITY 

The key microscopic notion behind the RFOT method- 
ology is that, regardless of the detailed interparticle po- 
tentials, local aperiodic arrangements in classical conden- 
sates become metastable below a certain temperature Ta 
(or above a certain density) [H, [26[ . Chemical detail and 
molecular structure affect the value of Ta, and the vis- 
cosity of the fluid. If the viscosity is high enough, one 
may cool the liquid so that it becomes locally trapped 
in metastable minima, while avoiding the nucleation of a 
periodic crystal, which would have been the lowest free 
energy state. Another system-dependent quantity is the 
size a of the elemental structural unit in the metastable 
liquid, or "bead" : the length a plays the role of the lat- 
tice spacing in the aperiodic structure, and is indeed quite 
analogous to the size of the unit cell in an oxide crystal, 
or it may correspond to the size of a rigid monomer or 
side chain in a polymer. The size a characterizes the 
range of the local chemical order that sets in during a 
crossover, at a temperature T cr , from collision dominated 
transport to activated dynamics [27| . The temperature 



F(N) 




(a) 



FIG. 1: (a): Typical nucleation profile of one aperi- 
odic lattice, within another, in deeply supercooled liquids. 
N = (47r/3)(r/a) 3 , AT* is the typical transition state size: 
(dF I dN) N t =0. £ is the volumetrically denned cooperative 
length: iV* = (£/a) 3 , where F(N*) = 0. (b): Cartoon of a 
structural rearrangement. The shown magnitude of £ corre- 
sponds to a temperature near T g on 1 hour scale. The two 
sets of circles - solid and dashed ones - denote two alternative 
structural states, dh — a/10, or the "Lindemann length", is 
the typical bead displacement during a transition. 



T cr is related to the meanfield temperature Ta but is al- 
ways smaller. The bead size a may be unambiguously 
determined from the fusion entropy of the corresponding 

271 ] . or else 



crystal, when the latter entropy is known 
can be computed from the fragility D using the univer- 
sal relationship between the latter and the heat cap acity 
jump at T g : D = 32/Ac p , as derived in RFOT @, Hf. 
Alternatively, if the configuration entropy can be reli- 
ably estimated, one may use the RFOT-derived relation 
for the configurational entropy per bead s c ~ .8fce 0]j 
which is somewhat sensitive to the barrier-softening ef- 
fects though [13 • 

Once locally metastable, the liquid may reconfigure 
but in an activated fashion, i.e. by nucleating a new ape- 
riodic structure within the present one. Such activated 
events occur, on average, once per typical a-relaxation 
time t, per region of size £. The nucleus grows in a 
sequence of individual, nearly-barrierless bead moves of 



length dh ~ a/ 10 [25|, [29( and time r„ 



I ps [21 



The overall sequence of elemental moves typically corre- 
sponds to the following activation profile, see FigQJa): 



(1) 



F(N) = 7 VN - Ts c N, 



where N is the size of a reconfigured region. The "sur- 
face term" jy/N is the mismatch penalty for creating 
one aperiodic structure within another. s c is the excess, 
"configurational" entropy of the liquid per bead, hence 
the entropic, bulk term (—Ts c N), which drives the tran- 
sition and reflects the multiplicity of possible aperiodic 
arrangements in a region of size N. The maximum of the 
profile: 



F* = m&x{F(N)} 
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(2) 



is achieved at N$ = N*/A, where F(N*) = 0, so that the 
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typical relaxation time is 

T -T . p F % /k B T = . DT K /(T-T K ) / o\ 

This formula works well at t > 1 nsec or so. The form on 
the r.h.s. is the Vogel-Fulcher law, derived in the RFOT. 

The end result of a cooperative, activated event is a 
reconfigured region of size £, where each of the N* = 
(£/ a ) 3 beads has moved the Lindemann length d^, or 
so, see FigQJb). Both £ and the nucleation critical size, 
r$ = C/4 1 / 3 , increase with lowering the temperature, 
roughly as r* cx £ cx 1/(T - T K f' 3 [J H i|. Here, 
Tk is the so called ideal glass transition temperature, 
where the excess liquid entropy s c , extrapolated below 
T g , would vanish [30J. At Tg, £ is still quite modest, 
only about six beads across [2|, [6| . Activated transport 
becomes dominant below the temperature T cr , such that 
r*(T CI ) = a, which corresponds, apparently universally, 
to r/r m i cro ~ 10 3 , or viscosities on the order of 10 Ps [27j |. 
At times shorter than 10 3 r m i C ro, one may then speak of 
a local aperiodic lattice on length scales of the cooper- 
ativity size £, since the slow structural reconfigurations 
have now time-scale separated from the vibrations [27l ] . 
Because of the local nature of structural relaxations, one 
speaks of dynamic heterogeneity, or a "mosaic" of coop- 
erative rearrangments 0] ■ The heterogeneity is two-fold: 
On the one hand, a local rearrangment implies that the 
surrounding structure is static during the transition, up 
to vibrations. On the other hand, because of the spatial 
and temporal variation in the local density of states (and 
hence variations in s c ), local reconfigurations are gener- 
ally subject to somewhat different barriers in different 
regions [2| (see Eq. © and also Section HV|) . 

Computation of the viscosity in such a dynamically 
heterogeneous environment may be done in two steps: 
First compute the viscosity in a medium with a homo- 
geneous relaxation time, call it r', and then average out 
with respect to the true distribution of the relaxation 
times. This procedure is valid in view of the equiva- 
lence of time and ensemble average. When the relax- 
ation rate is strictly spatially homogeneous, one may for- 
mally define a diffusion constant for an individual bead: 
D' = oj^/Gt 1 , since a bead moves the Lindemann length, 
once per time r', on average. Note that since a bead's 
movements, as embodied in dh and t', are dictated by 
its cage, this is an example of "slaved" motion, to bor- 
row Frauenfelder's adjective for conformational changes 
of a protein encased in a stiff solvent [3l|, H2] • One may 
associate, by detailed balance, a low-frequency drag co- 
efficient to that diffusion constant: (' = ksT/D'. Such 
dissipative response implies irreversible momentum ex- 
change between a chosen particle and its homogeneous (!) 
surrounding, hence a Stokes' viscosity rj — £'/(67ra/2), 
where a/2 is used for the radius of the region carved out 
in the liquid by a single bead. Averaging with respect 
to t' yields for the steady state viscosity of the actual 
heterogeneous liquid: 

V = -ZZjT ( r > . ( 4 ) 



where we have removed the prime at t, implying averag- 
ing with respect to the actual barrier distribution. The 
equation above can be rewritten as 

v =(a/d L r^f(T)^60^(r), (5) 

since the Lindemann ratio, dhja, ~ 0.1 has been argued 
to change at most by 10% between T m and T g [29| . 

Another instructive way to present Eq.(0J is to note 
that the Lindemann length is nearly equal to the typi- 
cal amplitude of high-frequency vibrations: — d V ib r , 
within 5% or so, see Fig. 3 of Ref.[29|. The vibrational 
amplitude is fixed by the equipartition theorem, since 
per bead: -fCx)a 3 (dvibr/a) 2 = &bT, where Koo is the high- 
frequency elastic constant of the aperiodic lattice. One 
gets, as a result, a Maxwell-type expression: 

V ~ (r) . (6) 

The last equation provides an easy way to see that the 
estimates in Eqs.(|4])-([6]) agree well with the experiment: 
Judging from the sound speed in glasses [H, [34|, the 
typical high frequency elastic modulus is about 10 9 — 10 10 
Pa, i.e. comparable but somewhat less than those of 
crystals. The range of relaxation times 10~ 12 — 10 4 sec, 
implies 10~ 3 — 10 13 Pa-sec for the viscosity, as is indeed 
observed. Alternatively, one may obtain these figures by 
substituting a typical a ~ 3A (see [13, [H| for specific 
estimates of bead sizes/densities). 

Finally, the exploited equivalence between the time 
and ensemble averages implies that crystallization has 
not begun during the experiment, of course. The latter 
possibility adds uncertainty into viscosity measurements, 
as the presence of crystallites would greatly broaden the 
dynamic range of local relaxations owing to relatively 
slow crystal nucleation events and the slow hydrody- 
namics near the crystallites. Similarly, long chain mo- 
tions in polymeric melts would also introduce additional 
long time scales into the problem. Our derivation does 
not apply to those situations. We note that optical 
transparency, which is often used as an indicator of no- 
crystallinity, does not ensure that crystallites - hundreds 
of nanometers across or smaller - are absent. Therefore 
a rigorous experimental study should, in the least, check 
whether performing viscoelastic measurements has en- 
hanced the crystallization of the sample. Ideally, X-ray 
diffraction should be monitored in the course of viscosity 
measurements. 



III. IONIC CONDUCTIVITY 

In any supercooled melt, whether regarded ionic or 
not, the beads carry an additional charge, relative to the 
corresponding crystal, because of the lack of crystalline 
symmetry. As a result, each structural reconfiguration is 
characterized by a transition-induced electric dipole, see 
FigH 
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construction of liquid states [f|. Repeating that argu- 
ment, but in the presence of electric field E, yields for 
the typical free energy profile for structural reconfigura- 
tion in steady state: 



F(N) = 1 VN~ t i T (N)E c -Ts c N, (10) 
where N is the size of the rearranged region and 



N/2 



(in) 



(11) 



FIG. 2: Shown on the left is a fragment of the mosaic of 
cooperatively reconfiguring regions in the supercooled liquid. 
Expanded portion shows how rotation of a bond leads to gen- 
erating an elemental dipole change during a transition, where 
the partial charges on the two beads are ±£g. 



The latter was estimated to be about a Debye or so, for 
most molecular substances [35| . This value comes about 
as we may break up the whole domain into N* jl pairs, 
where is pair has an elemental dipole fii, i = 1 . . . N* /2: 
(|jUj|) = (eq. Here q is the elementary charge: q = 10°e, 
and C < 1 characterizes the excess charge. This quantity 
£ is usually small reflecting small deviations from the 
crystalline symmetry, in the case of molecular crystals, or 
reflecting the weak interaction in Van der Waals systems. 
Alternatively, the overall density of charged/polar beads 
may be low. Ionic melts, by the very meaning of the 
term, are distinct from molecular/ Van der Waals systems 
in that nearly all beads are strongly charged, implying 
£ ~ 1. To be more specific, the conclusions of this article 
will be exemplified with an often studied mixture of 40% 
Ca(NO 3 ) 2 -60%KNO 3 ("CKN"), T g ~ 330K. 

During a transition, each dipole turns by an angle 
(di/a). The total transition dipole: 



N'/2 

i* = E WP 



(in) 



(7) 



scales as V N* because of the random orientation of the 
elemental dipoles (3o| : 



(IT 



CM[(£/ a ) 3 /2] 1/2 (dL/«). 



(8) 



When the dipole density is uniform, every transition re- 
sults in a local arrangement equally representative of the 
liquid structure. In other words, structural transitions 
do not modify the overall pattern of the immediate coor- 
dination shell. Transitions lead to a local ionic currents: 
i' = fi' T /T', per region of volume £ 3 . In the presence of 
an electric field, the net current density: 



3 = 



(9) 



is non-zero because the dipole moment at the transition 
state is correlated with the overall transition dipole mo- 
ment. The latter can be shown using Wolynes' library 



is the overall dipole change in that region. The subscript 
"c" in E c signifies that the latter is a cavity field. The 
field dependent term \/j,t(N)E c \ <C k B T is overwhelm- 
ingly smaller than the other terms in Eq. (Ti"(3| . As a result, 
the transition state dipole moment is not field-induced, 
but, again, is "slaved" to the lattice. More formally, one 
may use the argument from Ref. Q showing that the den- 
sity of structural states at the reconfiguration bottle-neck 
is of the order 1/T g , implying the field will not affect 
the specific sequence of elemental moves, but will affect 
the dynamics merely by shifting the energies along struc- 
turally dictated sequences of moves. Thus in the lowest 
order in E c , t- x {E c ) = t-\E c = 0)(1 + ^ T E c /k B T), 
yielding 



H T {H X T E C 



(12) 



where n\, = /j,t(N$) is the transition dipole moment at 
the zero-field transition state. The cavity field, see e.g. 
[3(J, is related to the external field E^- 00 ^ by 



E c ^) = E^)-^l 



(13) 



where ef,(o->) is the dielectric constant of the surrounding 
bulk. Since a steady current is implied in the derivation, 
(the imaginary part of) eb(u>) diverges at zero frequency, 
implying E c = E^-°°\3/2). One thus obtains for the ionic 
conductivity tensor: 



2k B re 



(14) 



Bearing in mind that A* = A*/4 = (£/a) 3 /4, and 
that the liquid is isotropic, on average: ((A/x)j(A/Lt)j) = 
Sij(Afi) 2 /3, one finally has: 



■'• '>><rknl\-r 



where 



A Mmol = 



,(/) 



(in) 
l 3 



(15) 



(16) 
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is the average elemental dipole change squared. Finally 
note that in covalently networked materials, where dipole 
assignment may be ambiguous, one may still estimate lo- 
cal dipole changes using the known piezoelectric proper- 
ties of the corresponding crystal, see [35| for details. 

The derivation above does not apply to systems where 
the dipole density is significatly non-uniform. For in- 
stance, glycerol has one polar, OH group per non-polar, 
aliphatic group, implying the liquid is non- homogeneous, 
dipole moment wise, on the a-relaxation time scale. In 
such systems, the premise that structural rearrangements 
result in equally representative configurations of the liq- 
uid does not hold. In the glycerol example, ionic conduc- 
tion would imply breaking OH or CH bonds. In CKN, 
on the other hand, the overall bond pattern, around any 
atom, does not change significantly during a transition, 
even though individual bonds distort by the Lindemann 
length, as mentioned. Eg . ([Toll thus places the absolute 
upper limit on the intrinsic ionic conductivity of a melt. 
By "intrinsic" we mean that the computed currents are 
always present in the fluid and result from the intrin- 
sic activated transport: Local bond pattern does not 
change significantly in the course of an individual ac- 
tivated event, but only in the course of many consecutive 
events, since during an individual event, the molecular 
displacements barely exceed typical vibrational displace- 
ments. Conversely, if a system displays a higher con- 
ductivity than prescribed by Eq. (|15p , one may conclude 
that the ion motion does not require structural reconfigu- 
ration. Here, the bond pattern actually changes, however 
these are not scaffold bonds of the aperiodic lattice com- 
prising the fluid (or glass). (More on this below.) 

To simplify comparison of Eq.(|15p with experiment, 
let us express the combination of the bead charge £q 
and size a, in Eq. (|15p . through the finite frequency di- 
electric response, a measurable quantity in principle (see 
below) . The latter is the response of a rearranging region 
in the absence of bulk current, i.e. with a fixed environ- 
ment, up to vibrations. It is convenient to choose such 
regions at volume £ 3 , so that each region has two struc- 
tural states available, within thermal reach from each 
other, separated by a barrier sampled from the actual 
barrier distribution in the liquid. If the two states, "1" 
and "2" are characterized by dipole moments fa and 
fa respectively, the expectation value of the dipole mo- 
ment of the region is /i = (fa + fa)/ 2 + A/j,(p2 ~ Pi)/2, 
where A/x = (fa — fa); p\ and p% = (1 — pi) are the 
probabities to occupy state 1 and 2 respectively. The 
relative population (j>2 — Pi) depends on the field via 
Sln(p2/pi) — AfiE c / ksT . At realistic field strengths, 
i.e. \AfiE c /kBT\ <C 1, one has for the field-induced shift 
of the relative population: 5(p2 — pi) — 2pip2(A[iE c ). 
Similarly to the preceding argument, ((Afi)i(Afi)j) = 
(iV*/2)<5ij(A/i mo i) 2 /3. Further, since we have frozen the 
structural transitions in the surrounding region, in es- 
timating the cavity field, one must use e(ui) with the 
a-relaxation contribution subtracted. This does not 
introduce much ambiguity because in most ionic sub- 



stances, the dielectric constant even at very high fre- 
quences is significantly larger than unity. In CKN, for 
instance, e'(oo) ~ 7 allowing us to write as before: 
E c ~ £K°°)(3/2). One thus obtains, in a standard fash- 
ion, for the frequency dependent dielectric constant in 
the absence of macroscopic current: 

ei„sM-eoo=Mp lP2 )-^f^/-^— V (17) 

Aa^KB-l \ 1 — IWT / 

where the label "ins" signifies the absence of dc conduc- 
tivity. 

In the presence of steady current, the full response per 
domain is the sum of the steady current from Eq. (|15|) 
and the ac current from Eq. (fP7|) . The addition of the 
dc contribution, ia/u, to the full dielectric response will 
increase the absolute value of e&(w). This means that 
Eq. (fTT|) should work even better. One thus gets for the 
full dielectric response of a conducting substance: 

e(u;)-e OQ =4n(p 1 p2)^^-(-^—)+ l -, (18) 

where a is the dc conductivity from Eq. (fT5)) . 

One needs to know the distribution of the transi- 
tion energies E to estimate the quantity (P1P2) = 
(l/4:COsh. 2 (E/2kBT)). Since £ is the smallest possible 
size of a rearranging unit, these rearrangments corre- 
spond to the elementary excitations in the system. We 
thus conclude, based on equipartition, that the typi- 
cal value of E is roughly fc^T, implying that (P1P2) is 
close to its maximum value of one quarter but is likely 
smaller by another factor of two or so. Assuming then, 
for the sake of argument that (P1P2) = 1/8, one gets 
A/1 2 

a j fc m y ~ 2(e' — eoo), within an order of magnitude. By 
Eq. (fTS]) , this implies a Maxwell- like relation between the 
dc conductivity and the real part of the dielectric re- 
sponse: 

*~(eo-eoo)(^). (19) 

with an important distinction, though, that here one av- 
erages the inverse relaxation time. The (e' — ) dif- 
ference in CKN, to be concrete, is about (20 — 7) = 13 
10, HZ!. This implies, by Eq.^ and CKN's T g ~ 330K 
37j ] . that (q ~ 3e, at a = 3A, a reasonable value for the 
bead charge. Note that a = 3A is consistent with CKN's 
Ac p ~ .12 cal/g K [38| and the mentioned D = 32/Ac p 
0. 

Now, substituting CKN's (e' — £00) into Eq. (fT5")) yields 
for the conductivity a ~ 3 • 1CP 10 (Ohm m) _1 sec. 

Naively replacing (t -1 ) with 1/ (r) would imply, at the 
glass transition, where (t) ~ 10 2 — 10 3 sec, a conductivity 
of the order 10~ 13 — 1CP 14 (Ohm cm) -1 , which is three 
to four orders of magnitude below the observed value 
[H H3- Note that the value (t -1 ) (t) = 10 3 - 10 4 is 
just the magnitude of decoupling observed in CKN near 
T g [ll], [k| , and is in fact expected for a fragile substance 
such as CKN is, as we will argue in the following. 
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IV. BARRIER DISTRIBUTION AND THE 
DECOUPLING 

Relaxation barriers in supercooled liquids are dis- 
tributed because the local density of states is non- 
uniform, leading to variations in the local value of the 
configurational entropy and hence the RFOT-derived 
barrier from Eq.([2]). In the simplest argument, the gaus- 
sian fluctuations of the entropy translate into gaussian 
fluctuations in the barrier, where the relative deviations 
of the two quantities, from the most probable value, are 
given by 



SF 



SF 



Ss c 

Sr 



1 



2-JD 



(20) 



where D is the liquid's fragility from Eq.Q. The quan- 
tity \j2\l~D varies between 0.05 and 0.25 or so, for known 
glassformers, the low and high limits corresponding to 
strong and fragile substances respectively. 

Xia and Wolynes (XW) further argued that the real 
barrier distribution should be cut-off at the most proba- 
ble value because a liquid region with relatively low den- 
sity of states is likely neighbors with a relatively fast 
region 3]. In addition we may recall that in the li- 
brary construction, the most likely liquid state is the one 
where the liquid is guaranteed to have an escape trajec- 
tory This means that the most probable barrier is 
also the maximum barrier. One may conclude then that 
the naive Gaussian distribution is adequate at small bar- 
riers, but significantly overestimates the probability of 
barriers larger than the typical barrier. Put another way, 
the trajectories corresponding to higher than most prob- 
able barrier in the naive Gaussian, all contribute to the 
F < F mp range. XW have implemented this notion by 
replacing the r.h.s. of the simplest Gaussian distribution 
by a delta-function centered at F mp [3J: 



Pi(F) = 



,-(l/_F-l) 2 /2SF 2 



2tt(SF) 2 F 2 



+ -S(F-1), (21) 



where F = F/F mp < 1 + , and we took advantage of 
the temperature-independence of the relative width in 
Eq. (|2H)l . This approximate form does not use adjustable 
parameters and quantitatively accounts for the correla- 
tion between the fragility and the stretching exponent 
(3 and the deviations from the Stokes-Einstein rela- 
tion. The distribution in Eq. (|2"Tj) is shown in FigJ3J The 
only difference of Eq. pTj) with the XW's form is that 
they used a purely gaussian form for F < 1, whereas 
we follow their own suggestion and employ the more ac- 
curate F cx l/s c (where s c is gaussianly distsributed 
of course). The accurate evaluation of the left wing of 
the distribution is imperative in estimating the average 
rate T~ i 1 cro e _F / fesT , because the latter is a rapidly vary- 
ing function of F. (&_bT is significantly less that SF at 
low temperatures.) Note that because of the rapid de- 
cay of the exponential at small F in Eq. (f2"T]) , accounting 




FIG. 3: The barrier distributions from Eqs.([2T} and ([23}. 
The S- function portion of the distribution from Eq. (|21[) is not 
shown. SF = 0.25. 



for the lowest order, quadratic fluctuations of entropy 
suffices. The quantity (r) that characterizes the 

apparent decoupling, computed with the XW's distribu- 
tion, is shown with the dashed line in FigHJ at T g , as a 
function of the relative distribution width SF. 

How robust is the prediction based on the simple func- 
tional form for the barrier distribution from Eq. (l2~T1) ? In 
spite of its quantitative successes, one may argue that 
the true barrier distribution should be a smoother func- 
tion, near F = 1. One way to see this is to computing 
e"(uj) from Eq.(JT7|» via the distribution in Eq.lpT]): The 
obtained curves are a sum of two peaks, one of which is 
broader, one the other is narrower than the experimental 
e"(u>)- The two peaks correspond to the half-Gaussian 
and the delta-function in Eq. (|2Tj) respectively. Let us see 
that knowing the precise form of the barrier distribution 
however is not essential in quantitative estimates of the 
decoupling so long as we account correctly for the overall 
width of the distribution and its decay at the low barrier 
side. 

It is straightforward to show that there exists a distri- 
bution that (a) satisfies these requirements without intro- 
ducing adjustable constants, (b) reproduces the experi- 
mental e(w) and does as well as the XW form for the f3 
vs. D correlation. As we have already discussed, the low 
barrier wing of the distribution in Eq. (|21[) is adequate. 
On the other hand, the high barrer wing should include 
the contributions from both sides of the original Gaussian 
peak, which are both of width SF/2. "Stacking" these 
two on top of each other, to the left of F mp , results in a 
distribution of width SF/4 (see aslo Appendix). Further, 
based on the known e{ui) data, the barrier distribution 
should be well approximated by an exponential, suggest- 
ing we use p(F) cx e F ^ SF ^ 4 ^ near F mp . Indeed, this im- 
plies p(r) cx T ( Ak BT/5F)-i _ frequencies not too close 
to the maximum of e"(w) and the rapid drop-off at small 
F, one has an approximate power law: 



i + H 5 



(22) 
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FIG. 4: The decoupling between the viscosity and the in- 
trinsic ionic conductivity, as a function of the relative barrier 
width 8F from Eq.([20j), at T 9 . The dashed and solid line 
pertain to the specific barrier distributions from Eq. (|21|l and 
(|23H respectively. 



FIG. 5: Different relaxation times derived from the bar- 
rier distribution in Eq, (|23|l . as functions of the most prob- 
able barrier (left) and the corresponding temperature (right). 
SF = 0.25, corresponding to fragility D = 4. 



We thus arrive at the following form: 



p(F) = l f 2 ~ - ~ Z 

\f 2 e F /( SF ^\ F e <F<l, 



(23) 



where F e and the normalization constants c\ and C2 are 
chosen so that the distribution is normalized, continuous, 
and its first derivative is continuous too. The distribution 
from Eq. ([2"3"|) is plotted in FigJ3J The decoupling strength 
(t) (V -1 ), computed for the composite distribution from 
Eq. (f2"3")l . is shown in Fig(U as a function of the relative 
distribution width SF, at T g . We therefore observe that 
in fragile liquids, the apparent time-scale separation may 
reach as much as four orders of magnitude near the glass 
transition - even though only one process is present! - 
because the inrinsic ionic conductivity is dominated by 
the fastest relaxing regions. 

Conversely, when the apparent decoupling exceeds the 
intrinsic value prescribed by Fig[4j we may conclude that 
ionic conduction does not in fact require structural relax- 
ation. This notion is of significance for the mechanisms 
of electrical conductance in glasses and will be discussed 
in detail in the Conclusions. 

One may also illustrate the effects of apparent decou- 
pling for a specific value of fragility, by plotting several 
varieties of relaxation times, as functions of the most 
probable barrier, or the corresponding temperature, see 
FigEl (Given the time scale at the glass transition, say 
DTk I (T g - T K ) = ln(10 16 ) ~ 37 (see Eq.©), there is 
a one-to-one correspondence between the fragility D and 
the T 9 /Tk ratio.) We observe that the average relax- 
ation time and the one derived from the inverse of the 
maximum position of e"(u>) are close, and are near the 
most probable value of the relaxation time. (The e"(u>) 
was computed using Eqs. fTg)) and (|2"3")l . see below.) The 
apparent conductivity relaxation time is strongly decou- 
pled, consistent with data of Howell at el. [ll| for CKN. 



Note that the value of fragility used in Figj5l D = 2, 
is probably smaller than in CKN. In addition, we have 
ignored here, for clarity, the effects of barrier softening 
[27j , that would require introducing a system-specific ad- 
justable constant Ta- The latter effects would change the 
slopes of the curves somewhat, without affecting their 
vertical separations. 

To test the predictions from FigsfJ] and O one needs 
to know the width of the barrier distribution for a- 
relaxation. As already mentioned, the gross features 
of this distribution have been predicted by the RFOT 
theory, and have lead to quantitative predictions of the 
correlation between the stretching exponent (3 and the 
fragility D, and the deviations from the Stokes-Einstein 
relation. The corresponding trends are as follows: more 
fragile liquids are predicted to have broader barrier distri- 
bution leading to a smaller value of [3, and vice versa for 
stronger substances [3j. A correlation with the fragility 
comes about by virtue of Eq. (|20p . Several d priori ways 
to determine and D have been employed, by experi- 
menters, that sometimes produce conflicting results. For 
example, the fragility extracted from r a will be consis- 
tently lower than that extracted from the mechanical 
relaxation time t s , because r CT < r s . The exponent 
from the stretched exponential is extracted from fits of 
various relaxation processes to a stretched exponential 
profile e~(*/ T )' 3 . Alternatively, one may choose to fit the 
Fourier transform of the stretch exponential, or the Cole- 
Davidson form, to the imaginary part of e(u) in insula- 
tors . These usually produce comparable results for 
the corresponding exponent (3, with a notable exception 
of ionic conductors, which happen to be the main focus of 
this paper. In ionically conducting systems, the dc com- 
ponent of the full dielectric response from Eq. (fT5|) largely 
"swamps" the ac part so that reliable determinations of 
the latter are complicated. The reader is reminded that 
dielectric measurements on ion melts are difficult because 
electrodes generally block ionic current. The effects of 
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FIG. 6: The top panel shows the imaginary part of the 
non-conductive part of the dielectric response, e(u>), Eq, (|17[l . 
for four values of the most probable a-relaxation barrier: 
Fmp/fcsT = ln(r/r m i cro ), as indicated on the graph. (SF = 
0.25, 



FIG. 7: The widths of the peaks in the imaginary parts of 
the dielectric constant e(w) and modulus M{ui), as functions 
of temperature (c.f. Fig[6j. SF = 0.25. 



We have used CKN's values 
For the sake of argument, wc 



distribution from Eq 
for 6q and e^, as before 



Tmir-rn = 1 PS.) These barrier values were chosen be- j p noc j- j- a n /i + t t?- m 

micro 7 1n3 . use oi* = 0.25, corresponding to p ~ 0.4 at i ff . lnlig|6J 



cause for T/r m i cro > e 7 — lO 3 , the RFOT approach is quan 
titatively accurate, while r/r ro i cro — e 37 ~ 10 16 is at the 
upper limit of the dynamical range routinely accessible in the 
lab. The bottom panel show the corresponding modulus, from 
Eq. (|24[l . that includs the dc-part, due to the intrinsic ionic 
current. The same four values of the barrier are used. We 
have used the CKN's values for eo and eoo, and assumed that 



A 1.1 



mol 



see the discussion preceeding Eq. (fl9 



build-up charge are often treated phenomenologically, by 
means of equivalent circuits [3?], Eo] • Given these compli- 
cations, many have chosen to p lot the reciprocal of e(u>), 
i.e. the dielectric modulus [ill. [40|: 



M(cj) ee l/e(w). 



(24) 



M(ui) is well behaved and even shows a peak in the 
imaginary component, similarly to e(u>) of a near insu- 
lator. In the absence of an a priori microscopic picture 
and by analogy with e(w), one might interpret this peak 
as as the response of the electric field E to the dielec- 
tric displacement D. This in fact would be appropriate 
in a layered dielectric [Till . See also the discussions in 
Refs.fll^HIllHEnil. Yet the resulting values of 
the most probable relaxation time and the stretching ex- 
ponent deviate from those obtained with other methods 
[i~3[ fig. l37|. In fact, the modulus-derived /3m increases, 
while the width of the M"(u) peak decreases with lower- 
ing the temperature, in conflict with the general trends 
for poor conductors, and the conclusions of the RFOT 
theory. 

The RFOT theory and the present results allow one 
to address these difficulties, to which we devote the rest 
of this Section. One first notes that structural reconfig- 
urations are compact, and so the layered-dielectric view 
of supercoold melts is not microscopically justified. We 
next plot, in the top panel of Fig[51 the non-conductive 
e[' ns (uj), from Eq. fTT)) averaged with respect to the barrier 



bottom, we show the imaginary part M"(u>), of the full 
modulus. Clearly the two functions exhibit qualitatively 
different behaviors. Note that the effect of the dc com- 
ponent on the apparent relaxation profile has been dis- 
cussed previously [13, 01 , including the possibility of a 
double peak [24j . The latter has been observed by Funkc 
at e l. 1 411) . but has not been reproduced by Pimenov at 
el. [33] • Nevertheless, the dielectric modulus obtained 
here is qualitatively consistent with CKN's data from 
Ref . [37| . Finally note that for smaller dc-conductivities, 
the modulus data would become more similar to e"(cj). 

We conclude from the above analysis that if one were 
to use the modulus data to extract the characteristics 
of the barrier distribution, one must measure first the 
dc-current, add it to the ei ns (uj) from a microscopic the- 
ory, and then compare the result to the measured M(uj) 
data. But again, because of the large contribution of 
the dc component, the corresponding fits would not dis- 
criminate well between different forms of ei ns (w). On the 
other hand, treating the electric field as a response to the 
displacement may lead to erroneous conclusions on the 
temperature dependence of the barrier width. In fact, 
the barrier widths derived from e(n S ( w ) or M"(u) show 
the opposite trends, as we have seen already in Figj6] 
One may further quantify this observation: In the ab- 
sence of a microscopic theory, one often characterizes the 
width of the 6-^(0;) peak by a stretching exponent /?, 
as derived e.g. from Davidson-Cole fits. The distribu- 
tion from Eq. (|23[> indeed gives rise a power law behav- 
ior, consistent with Eq. l|2"2"|) , see Appendix. In contrast, 
the corresponding M"(lu) curves do not exhibit a similar 
power-law behavior. I have chosen to illustrate the oppo- 
site temperature trends in the widths of e-^ s [ln(o;)] and 
M"[ln(w)] peaks, by measuring the latter widths at one- 
third-height and plotting them as functions of temper- 
ature, see Fig[7] Similar opposite trends, too, would be 
observed for the corresponding apparent barrier widths or 
effective /3's. Clearly, interpreting the dielectric modulus 
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of an ionic conductor as a response function may lead to a 
significant underestimation of the actual barrier width at 
low temperatures, and qualitatively incorrect conclusions 
on the temperature dependence of the width. 



V. CONCLUSIONS 

We have computed, from the first principles, the vis- 
cosity and the intrinsic ionic conductivity of supercooled 
liquids. The viscosity is determined by four microscop- 
ically defined quantities: the length scale of the local 
chemical order that sets in at temperature X^, where liq- 
uid dynamics become activated; the Lindemann length, 
characterizing molecular displacements at the mechanical 
stability edge; the temperature; and the average relax- 
ation time r of the activated reconfigurations that dom- 
inate the liquid dynamics below T cr . The extraordinar- 
ily long r range is what gives rise to the high viscos- 
ity of the liquid when it approaches the glass transition. 
When the local chemically stable units (or "beads") are 
charged, the fluid will also exhibit an ionic conductiv- 
ity, which we have called the "intrinsic" conductivity, to 
constrast it with electric conduction via delocalized elec- 
tronic carriers or via mobile ions that are not bonded 
to the metastable aperiodic lattice forming the super- 
cooled liquid. Computing the conductivity requires an 
additional microscopic characteristic, the electric charge 
on a "bead" . Fortunately, this additional parameter may 
be deduced from the ac dielectric response, which we have 
also estimated. Perhaps the main finding of this work is 
that in contrast with the viscosity, the ionic conductivity 
is dominated by the fastest relaxing regions in the liquid, 
as reflected in Eq. (fT5l) . 

We have discussed ways to test the above predictions, 
the most important aspect of which is the large separa- 
tion, or "decoupling" , between the apparent time scales, 
suggested previously by viscosity and ionic conductiv- 
ity data on purely phenomenological grounds. We have 
shown that such apparent time-scale separation is indeed 
expected because of the very broad barrier distribution 
for a-relaxation, derived earlier in the Random First Or- 
der Transition (RFOT) methodology. The decoupling 
thus stems essentially from the same cause as the vio- 
lation of Stokcs-Einstein relation in supercooled liquids 
4]. Now, we have seen that the value of the decoupling 
is not very sensitive to the precise form of the barrier dis- 
tribution so long as one acounts for the RFOT-derived 
gross characteristics of this. We have thus quantified the 
degree of "decoupling" : The intrinsic ionic conductivity 
was argued to decouple at most by four orders of magni- 
tude from the low-frequency momentum transport. Con- 
versely, any conductivity exceeding this limit must be 
due to other charge carriers that do not disturb the liquid 
structure beyond typical vibrational displacements. In- 
deed, suppose the apparent decoupling exceeds the value 
prescribed by the width of the barrier distribution. This 
means that there will be ions that travel a distance ex- 
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FIG. 8: The four thick lines are the same as in Fig(BJ but 
in the double-log scale. The thin lines are the Davidson-Cole 
forms, following from a simple approximation, see Appendix. 
The dash-dotted line illustrates how the stretching exponent 
/3 was extracted from the curves. 

ceeding the Lindemann length in a time it takes the local 
environment to relax. Therefore, local relaxation is not a 
necessary condition for a non-zero current of these ions. 
Some interaction with relaxation may still be present, 
however at large enough decouplings, we may say that 
the ion (or any other carrier) interacts with the liquid 
as if the latter were a perfectly stable, albeit disordered 
lattice. In such cases, one may think of the ionic current 
in superionic conductors in terms of regular, not slaved 
diffusion. In regular diffusion, the total travel time is 
dominated by the slowest step, in contrast to Eq. (fT5|) . 

The intrinsic difficulty in experimental assessment of 
the barrier distribution in moderately conductive melts 
is that the dc current dominates the overall dielectric re- 
sponse. This gives rise to ambiguities as to what the 
actual width of the barrier distribution is, since mechan- 
ical relaxation and dielectric modulus data disagree. We 
have shown that this is expected, and argued that the 
mechanical relaxation offers the preferred method of es- 
timating the actual barrier distribution. 
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Appendix 

Let us see that the distribution in Eq. (p?5|) is qualita- 
tively consistent with experimental ei ns (uj) and the em- 
pirical correlation of (3 and D. For this, we replot the 
top panel of FigHJ in the double- log format, in FigJS] 
We note the general adequacy of the barrier distribution 
from Eq. (|23p : Similarly to the experimental e(u>) in poor 
conductors, the resulting high-frequency wing is signif- 
icantly broader than the low-frequency one. Note that 
the actual data would also often display an additional 
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FIG. 9: Two approximate relations of the stretching exponent 
P to the fragility D from the Vogel-Fulcher form, from Eq. (|25p 
and as derived from the slopes of the high-frequency wing of 
the peaks, such as in Fig(S] 



fied by the dash-dotted line in Fig[8l The dependence 
of the thus obtained exponent (3 on the fragility D, at a 
fixed F mp /kBT = 37, is shown by the dashed-dotted line 
in FigEl This (3 is, again, qualitatively consistent with 
experiment. Greater accuracy should not be expected 
here, as we have not treated the higher- frequency range 
associated with /3-relaxation, which would affect the ex- 
perimentally determined stretching exponents. 

In addition, we verify that the informal argument in 
the main text that the width of the barrier distribution 
should be about 6F/A, at the half- height or so. Indeed, 
for a gaussian barrier distribution with width SF/4 = 
F/8^/D implies the following approximate expression for 
the stretching exponent (3 at T g (c.f. Eq.(9) of Ref.0]): 



1 + 



( F/k B T 



-1/2 



(25) 



high-frequency wing, which is ascribed to the secondary, 
/3-processes, also called Johari- Goldstein relaxation [421 ]. 
(See [HI for a review). The present results suggest that 
/3-relaxation does not contribute to the intrinsic ionic 
conductivity. At any rate, the derived e(co) show sev- 
eral decades of nearly power-law decay, allowing one to 
extract the corresponding exponent: e(w) oc uj^^ . The 
effective /?'s were deduced from the slopes of the curves 
at the points of maximum second derivative, as exempli- 



shown as the solid line in FigJSJ This expression is in 
very good agreement with experiment, see Fig. 2 from 
Ref.Q. (At T g on scale r/r micro = 10 16 , F/k B T ~ 37.) 
Note also Eg. ([23)1 is consistent with Eq.(|2"0l. assuming 
the Davidson-Cole [II] and William- Watts [45| stretch- 
ing exponents f3 are close [3j|. That the latter is the 
case indeed I demonstrate by graphing the Davidson-Cole 
(DC) form (e DO (w) - e^) = (e - eoo )(l - iujr)^ 
with r = T micro e Fm P /fcBT and (3 from Eq.fllS]). These are 
shown in FigJS] as thin dashed lines. 
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